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Abstract 

The Fermi bubbles are gigantic gamma-ray structures in our Galaxy. The physical origin of the bubbles is still 
under debate. The leading scenarios can be divided into two categories. One is the nuclear star forming activity 
similar to extragalactic starburst galaxies and the other is the past active galactic nucleus (AGN) like activity of 
the Galactic center supermassive black hole. In this paper, we propose that metal abundance measurements will 
provide an important clue to probe their origin. Based on a simple spherically symmetric bubble model, we find 
that the generated metallicity and abundance pattern of the bubbles’ gas strongly depend on assumed star formation 
or AGN activities. Star formation scenarios predict higher metallicities and abundance ratios of [O/Fe] and [Ne/Fe] 
than AGN scenarios do because of supernovae ejecta. Furthermore, the resultant abundance depends on the gamma- 
ray emission process because different mass injection histories are required for the different gamma-ray emission 
processes due to the acceleration and cooling time scales of non-thermal particles. Future X-ray missions such as 
ASTRO-H and Athena will give a clue to probe the origin of the bubbles through abundance measurements with their 
high energy resolution instruments. 
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1. Introduction 

The Fermi bubbles are gigantic gamma-ray structures ex¬ 
tending ~ 50° north and south of the Galactic center (GC) 
with a longitudinal with ^ 40° (Dobler et al., 2010; Su et al., 
2010; Ackermann et al., 2014). Structures roughly coincident 
with the gamma-ray bubbles are known in X-rays (Snowden 
et al., 1997; Bland-Hawthorn & Cohen, 2003), microwave 
(Finkbeiner, 2004; Dobler & Finkbeiner, 2008; Ade et al., 
2013), and polarized radio (Carretti et al., 2013). Past activities 
of our Galaxy is believed to generate these structures. At lower 
latitudes |6| ^ 20°, an additional gamma-ray emission compo¬ 
nent is reported in the bubbles (Hooper & Slatyer, 2013). This 
component would originate in millisecond pulsars or annihila¬ 
tion of dark matter particles rather than past activities of our 
Galaxy (Hooper & Slatyer, 2013), although the latest analysis 
of the bubbles does not find that component because of a large 
systematic uncertainty (Ackermann et al., 2014). 

Gamma-ray emission of the bubbles is thought to be sup¬ 
plied by leptonic or hadronic processes, namely the inverse- 
Compton scattering of interstellar radiation field and the cos¬ 
mic microwave background by electrons (e.g. Cheng et al., 
2011; Mertsch & Sarkar, 2011; Lacki, 2014) or the hadronu- 
clear process of protons (and ions) colliding with ambient gas 
in the bubbles (e.g. Crocker & Aharonian, 2011; Thoudam, 
2013; Fujita et al., 2013). Both models can explain the mi¬ 
crowave and gamma-ray data, although additional primary 
electrons or reacceleration of secondary leptons may be re¬ 
quired in the hadronic scenario (Fujita et al., 2014; Ackermann 


et al., 2014). 

In either case, the huge energy content of the bubbles an or¬ 
der of ergs (Su et al., 2010; Ackermann et al., 2014) 

should be explained as well. A fundamental question on the 
bubbles is what powers the bubbles. Theoretically, two scenar¬ 
ios are proposed as the origin of the bubbles. Those are nu¬ 
clear star-formation activity (e.g. Crocker & Aharonian, 2011; 
Carretti et al., 2013; Lacki, 2014) and past active galactic nu¬ 
cleus (AGN) activities of Sgr A* (e.g. Cheng et al., 2011; 
Zubovas et al., 2011; Guo & Mathews, 2012; Mou et al., 2014; 
Yang et al., 2013). Although a jet-like structure in the bub¬ 
bles was previously reported (Su & Finkbeiner, 2012), which 
supported the Sgr A* jet scenario, that structure was not con¬ 
firmed in the latest analysis (Ackermann et al., 2014). Carretti 
et al. (2013) has argued that the nuclear star formation activ¬ 
ity scenario is favored based on the polarization measurement. 
However, the measured polarization features have been argued 
to be also reproduced by the AGN jet scenario (Yang et al., 
2013). Other probes are necessary to investigate the origin of 
the bubbles. 

Kataoka et al. (2013), Tahara et al. (2015), and Kataoka et al. 
(2015) have recently carried out X-ray observations of the bub¬ 
bles using the X-ray Imaging Spectrometers (XIS) onboard the 
Suzaku X-ray satellite. The observed diffuse X-ray emission 
shows the existence of kT ~ 0.3 keV thermal plasma which 
is slightly hotter than the surrounding Galactic Halo (GH) gas. 
Tahara et al. (2015) have further found the possible existence 
of 0.7 keV plasma is indicated in the northern cap region which 
is seen in the all sky map of the Monitor of All-sky X-ray Image 


2 


[Vol., 


(MAXI). They found the expansion velocity of the bubbles as 
^ 300 km lower than most of previously proposed mod¬ 
els (e.g. Cheng et ah, 2011; Zubovas et al., 2011; Mertsch & 
Sarkar, 2011; Guo & Mathews, 2012; Lacki, 2014). This ve¬ 
locity is supported by the measurement of the X-ray absorption 
line toward 3C 273 whose sightline passes through the neigh¬ 
borhood of the bubbles (Fang & Jiang, 2014). Moreover, Fox 
et al. (2014) reported two high-velocity metal absorption com¬ 
ponents at -235 and h- 235 km/s using the spectrum of a quasar 
whose sightline passes through the bubbles. 

In this paper, we propose X-ray abundance measurements in 
the bubbles will provide a unique key to identify their origin 
because the distributed elemental abundances depend on yields 
of ejecta and mass loading factor of ambient gas. The region 
of the bubbles was initially filled with the low metal GH gas. 
In the star forming activity scenarios, the bubbles are polluted 
by the elements produced by supernovae (SNe) whose abun¬ 
dances are different from that in the interstellar medium (ISM). 
On the other hand, in the AGN wind scenario, the abundance 
of the wind would be the same as the ambient ISM which ac¬ 
cretes onto the Sgr A*. The resultant abundance distribution in 
the bubbles is expected to be different between the AGN wind 
and star forming scenarios. We also argue prospect for future 
X-ray observations. We adopt solar abundances reported in 
Asplund et al. (2009). Thus, the solar metallicity is set to be 
Zq ~ 0.0134 rather than classical value of Zq ~ 0.02 (Anders 
& Grevesse, 1989). 

In this Letter, we do not consider the the AGN jet scenario. 
The interior of the bubbles formed by jets would be polluted 
by metals of the jet itself. The Kelvin-Helmholtz instability 
is expected to be suppressed even with low level of viscosity 
(Guo et al., 2012). As the jets push the GH gas away, the metal 
mixing with the GH gas would not efficiently occur in the jet- 
induced bubbles. However, the jet composition is highly un¬ 
certain. Although pure pair jet models are excluded for blazars 
(Sikora & Madejski, 2000) and pairs may not survive the an¬ 
nihilation in the inner, compact and dense regions (Celotti & 
Ghisellini, 2008), there is still room for pairs in the jet, based 
on the energetics arguments (Sikora et al., 2005). Moreover, 
iron emission lines are observed in the jet of the Galactic mi¬ 
croquasar SS 433 (Migliari et al., 2002). 

2. Metal Enrichment in the Fermi Bubbles 

To consider the metal enrichment in the bubbles, first we 
consider the metallicity in the outflow. We follow the descrip¬ 
tions in Strickland & Heckman (2009), which discussed the 
outflow in the nearby starburst galaxy M 82. The net mass 
outflow rate from the GC is described as Mout = -^ejecta + 
Mism = PMejecta., where Mejecta N the ejected mass outflow 
rate from the origin to the bubbles, Mism is the loaded ISM 
mass rate, and /3 is the mass loading factor. If /3 = 1, no ISM 
gas is loaded. Given the star formation rate (SFR) and ini¬ 
tial mass function (IMF), the ejected mass outflow rate is esti¬ 
mated. Then, by comparing with the required total mass out¬ 
flow rate for the formation of the bubbles, the mass-loading 
factor is determined. For the nearby starburst galaxy M 82, the 
mass loading factor is in the range of 1.5 < /3 < 2.5 (Strickland 


& Heckman, 2009). 

The elemental abundance of an element i in the out¬ 

flow, i.e. the elemental mass fraction against the total baryons 
in the outflow gas, is 


( 1 ) 


where ejecta and ism is the abundance of the element in 
the ejecta and in the ambient ISM, respectively. Hereinafter, 
we assume X^ ism = (see e.g. Uchiyama et al., 2013; 
Nakashima et al., 2013, and references therein). 

Now we are interested in the abundance distribution in the 
bubbles. For the sake of simplicity, we assume a spherically 
symmetric bubble model. And, we also simply assume the GH 
gas had a distribution of pcH cx before the bubbles formed. 
Although the GH gas distribution has been under debate (see 
e.g. Yao et al., 2009; Miller & Bregman, 2013; Sakai et al., 
2014, for details), recent measurements by XMM-Newton sug¬ 
gest pgh cx at r ^ 0.35 kpc (Miller & Bregman, 2013). 
The adiabatic index of the gas is set to be 5/3. We adopt the 
self-similar solution for the hydrodynamical evolution of the 
gas (see e.g. Mihalas &. Mihalas, 1984; Ostriker &. McKee, 
1988). Depending on the material injection history, the result¬ 
ing matter distribution differs (see e.g. Fig. 2 in Fujita et al., 
2013). Instantaneous injection leads the gas mixing between 
outflow and the GH gas inside of the shock radius (i?sh), while 
the continuous injection leads a compressed GH gas between 
the shock and the contact discontinuity at i?cd = 0.84i?sh and 
only outflow gas exists behind Red- The metal abundance in 
the bubbles will be given as follows. Instantaneous injection 
case gives 

{ ^i,out Afouttwind-r^i.GHMGH (r < R i 

Afouttwind-|-*fGH ' ~ (2) 

Xi,GH {r > Rsh), 


while continuous injection case gives 


/ Xj out (r* ^ Red) 
\ Xi,GH {r > Red), 


(3) 


where Xi fb is the abundance of an element i in the bubbles, 
iwind is the time scale where the wind is active, X^ gh is the 
abundance of the element in the GH, and Mgh is the swept-up 
GH gas mass. The latter case is analogous to that in the wind of 
the starburst galaxies (e.g. Strickland & Heckman, 2009). We 
set Xi gh = 0.45Xi 0 (Miller & Bregman, 2014)*, although 
the GH gas metallicity is still uncertain (see also Sakai et al., 
2014, claiming solar metallicity). 

From the Suzaku observations, the shock radius is indicated 
at around 10 kpc from the GC (Kataoka et al., 2013). The 
swept-up halo gas mass is estimated as ~ 1.2 x 10 ®Mq us¬ 
ing the spherical /3 model (Miller & Bregman, 2013), while 
we assume the gas distribution follows r~^ for the simplic¬ 
ity. Once the abundances of the ejecta, the mass outflow rate, 
the timescale of wind activity and the mass loading factor are 


* We renomialize the reported value based on Anders & Grevesse (1989) to 
the latest solar abundance based on Asplund et al. (2009). 
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Table 1. Model Parameters for Metal-Enriched Outflows 


Origin 

Star formation 

AGN wind 


AGN wind 

Emission 

Leptonic 

Hadronic 

Leptonic 


Hadronic 

Reference 

Lacki (2014) 

Crocker et al. (2014) 

Mou et al. (2014) 

Zubovas et al. (2011) 

SER [Me/yr] 

0.1 

0.1 

- 

- 


IME model 

Salpeter (1955) 

Kroupa (2001) 

- 

- 


IME ranges 

0.1-100 Mq 

0.08-150 Mq 

- 

- 


Mont [ATg/yr] 

0.02 

0.1 

0.02 

0.08“ 


P 

2.0 

6.3 

_b 

_b 


ZfbIZq 

5.3= 

2 .2= 

1 .0= 

“045^ 


Afpe.Fs/ATpe,© 

2.3= 

1.3= 

1 .0= 

0.45'^ 


[O/Ee] 

0.49= 

0.30= 

0 .0= 

0 .0'^ 


[Ne/Ee] 

0.58= 

0.38= 

-FTTT-7^ 

0 .0= 

0 .0'^ 



b. 


i does not affect results assuming X^^ejecta = (see the details in the text). 


Expected values behind the contact discontinuity, Red- At larger radii, it will be the value of the GH gas. 
Expected values in the bubbles elsewhere. 


given, we can calculate the abundance distribution of the bub¬ 
bles from Eqs. 1, 2, and 3. 

In the nuclear star formation scenario, stars distribute ele¬ 
ments through SNe and stellar winds (SWs). We assume all 
stars have the solar abundances, since we assume ism = 
Xi^Q. In this paper, we neglect the yields of SWs, which may 
be crucial for light elements. We do not discuss the H-burning 
products below. The contribution of SWs to yields of heavier 
elements is expected to be small for stars having solar abun¬ 
dances even taking into account rotation (e.g. Hirschi et al., 
2005). Nomoto et al. (2006) provide the yields from various 
mass core-collapse SNe and hypernovae (HNe) whose explo¬ 
sion energy is > 10®^ ergs (Nomoto et al., 2006). The stars 
having mass of ^25-140 Mq in the main-sequence stage col¬ 
lapse to form a black hole. If the black hole has little angular 
momentum, little mass ejected. However, if the black hole ro¬ 
tates, the black hole eject matter through jet and it would be 
observed as a HN (Nomoto et al., 2013). 

We estimate the the SN ejecta abundances as follows 
(Nomoto et al., 2006). Given the IME (j){M)dM, the IME- 
integrated yields normalized by the total mass of ejected ma¬ 
terials are as follows (Nomoto et al., 2006; Tominaga et al., 
IQQlf: 

^ ImZT 

ejecta— 

C”“(Mej.SN[M] + Mej.sw[M])<^(M)dM 

where Xi ejecta is an integrated mass fraction of an element 
is mass fraction of i linearly interpolated between 
nearest models of Nomoto et al. (2006) as a function of an 
ejected mass, Mej^sN is an ejected mass by a SN, Mej^sw is 
an ejected mass by SWs, and M is the mass of a main se¬ 
quence star. Mijiin and Mmax is the minimum and maximum 
mass of stars, respectively. Eollowing Nomoto et al. (2006), 
We assume M < IOMq and M > 50Mq stars do not yield any 
materials, i.e. Mej^sN(Af < IOMq) = Mej 3 N(> 50Mq) = 0. 

^ In Nomoto et al. (2006), the IMF-integrated yields are normalized by the to- 
tal amount of gases forming stars. Since we are interested in the abundance 
in the ejecta now, we adopt the Eq. 4 in this Letter. 


We assumed the fraction of HNe to whole SNe ehn = 0 for 
M < 20Mq and chn = 0.5 for M > 20Mq (Kobayashi et al., 
2006; Nomoto et al., 2006). 

Eollowing Eq. 4, Xpe,ejecta for the Salpeter IME with the 
mass range of 0.1-100 Mq is 4.0 Xpe,©- In the nearby star- 
burst galaxy M 82, its outflow is predicted to have Xpe,ejecta ~ 
5A'Fe,Q (see e.g. Strickland &. Heckman, 2009), although the 
assumed IME and yields are different. 

In the case of past AGN-like activities of Sgr A*, the sit¬ 
uation is different. The ejecta abundances reflect the accre¬ 
tion disk abundances which are the same as the ISM abun¬ 
dances. Thus, we set Xi^ejecta = ATijsM = Xi^Q in the AGN 
disk wind scenarios. Eq. 1 implies that the yield of the out¬ 
flow is Wi out = at, Q. The mass loading factor does not affect 
results in the AGN disk wind scenarios. 

In this paper, we consider the leptonic star formation (SE) 
scenario (e.g. Lacki, 2014), the hadronic SE scenario (e.g. 
Crocker & Aharonian, 2011), the leptonic AGN wind (AW) 
scenario (e.g. Mou et al., 2014), and the hadronic AW scenario 
(e.g. Zubovas et al., 2011). The model parameters are summa¬ 
rized in Table. 1. As described below, we adopt the continuous 
injection case for the first three scenarios, while we adopt the 
instantaneous injection case for the hadronic AW scenario. 

Eor the leptonic SE scenario, we adopt the fiducial model 
parameters in Lacki (2014). They take the Salpeter initial 
mass function (Salpeter, 1955) ranging 0.1-100 Mq with the 
continuous SER of 0.1 Mq yr“^. The mass outflow rate is 
0.02 Mq yr“^ with j5 of 2.0. 

Eor the hadronic SE scenario, we adopt Crocker et al. (2014) 
where they adopt the Kroupa initial mass function (Kroupa, 
2001) ranging 0.08-150 Mq with the continuous SER of 
0.1 Mq yr“^ (Crocker, 2012). The mass outflow rate is set to 
be 0.1 Mq yr“^. The mass-loading factor is estimated as fol¬ 
lows. Given the SER and IME, the SNh-SW ejected mass out¬ 
flow rate is 0.016Mq yr“^. Then, /? = Mwind/Mejecta — 6.3 
assuming all the ejecta materials are injected into the bubbles 
(Crocker, 2012). 

Eor the leptonic AW model, we adopt the run A of Mou et al. 
(2014). They assume a radiative inefficiency accretion flow, but 
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2 x 10 ^ times higher accretion rate than present value motivated 
by Totani (2006) whose model can nicely explain various as¬ 
pects of the GC observables by past Sgr A* activity (see Totani, 
2006, for details). The accretion disk wind has the continuous 
mass outflow for 12.3 Myr. 

For the hadronic AW model, we adopt Zubovas et al. (2011) 
which assume an Eddington accretion wind but blowing only 
for fwind ~ 5 X 10^ yr at ~ 6 Myr ago. The mass outflow 
rate from the GC region is terminated in other epochs. Since 
the mass injection occurs for short time scale comparing to the 
age of the bubble, the hadronic AW model can be regarded 
as the instantaneous injection. As described in Zubovas et al. 
( 2011 ), the mass outflow rate is ^ 8 x lO“^M 0 yr“^ during the 
Eddington phase. 

3. Results 

The expected metallicity, iron abundance, and abundance ra¬ 
tios at a given radius are summarized in Table. 1. We note that 
the observed values are integrated values on the line of sight as 
a function of the Galactic longitude and latitude. The metallic¬ 
ity in the bubbles will be 5.3 Zq, 2.2 Zq, and Zq at r < i?cd 
for the leptonic SE scenario, the hadronic SE scenario, and the 
leptonic AW scenario, respectively. At r > R^d, it will be the 
GH gas metallicity. Therefore, as given in Eq. 3, the metal¬ 
licity in the bubbles would have a clear jump at the contact 
discontinuity at r ^ 8 kpc from the GC for the continuous in¬ 
jection cases. Because of the difference of the mass loading 
factor, the hadronic SE scenario predict lower metallicity than 
the leptonic SE scenario does. Since we assumed that the AGN 
disk wind and the loaded ISM have the solar abundance, the 
expected metallicity becomes Zq. Eor the hadronic AW sce¬ 
nario, it will be kept at the GH gas metallicity level, 0.45 ^ 0 , 
at elsewhere. Although there is a small metallicity jump at the 
shock radius, that will be a factor of < 0.5 % jump. This is 
because the injected gas amount ^ 4.0 x IO^Mq is relatively 
smaller than the swept-up GH gas mass 1.2 X 10®Mq. 

It is hard to distinguish models with current X-ray data 
through metallicities, since Suzaku data have huge uncertain¬ 
ties in deriving metallicities due to low photon statistics and its 
energy resolution. Eurther X-ray observations are required to 
unveil the origin of the bubbles through the abundance mea¬ 
surements. Interestingly, future missions such as ASTRO-H 
(Takahashi et al., 2012) and Athena (Nandra et al., 2013) will 
have high energy-resolution spectrometers, which may enable 
us to study abundance ratios. Once elemental line emissions 
are clearly measured, we can reliably determine the metallic¬ 
ities and abundances in the bubbles. To compare with future 
data, we also evaluate the iron abundance and the abundance 
ratios which are the logarithm of the ratio of abundances com¬ 
pared to the solar abundance ratio. The iron abundance in the 
bubbles behind the contact discontinuity will be 2.3 Xpe,©, 
1.3 Xpe,©, Xpe,© for the leptonic SE scenario, the hadronic SE 
scenario, and the leptonic AW scenario, respectively. The iron 
abundance in the hadronic AW scenario will be 0.45 Xpe,© at 
elsewhere. The abundance ratio of [O/Ee] behind the contact 
discontinuity will be 0.49, 0.30, and 0 for the leptonic SE sce¬ 
nario, the hadronic SE scenario, and the leptonic AW scenario, 
while it will be 0 for all models at larger radii. It will also 


be zero at elsewhere for the hadronic AW scenario. We note 
that the solar abundance ratio corresponds to zero. Thus, AW 
scenarios give the value of zero. [Ne/Ee] also give the similar 
results as in [O/Ee], but [Ne/Ee] will be 0.58 and 0.38 for the 
leptonic SE scenario and the hadronic SE scenario behind the 
contact discontinuity. 

We also perform spectral simulations for ASTRO-H^. Eigure. 
1 shows the simulated spectrum of the bubbles with 200 ks ex¬ 
posure for the Soft X-ray spectrometer (SXS) onboard AST/JG- 
H. Three components are included. Those are the bubbles, the 
local hot bubble, and the cosmic X-ray background following 
Kataoka et al. (2013); Tahara et al. (2015). Since the fore¬ 
ground GH gas component is not observed in the bubbles’ re¬ 
gion (Kataoka et al., 2013; Tahara et al., 2015; Kataoka et al., 
2015), the GH gas component is not included here. We assume 
the same spectral parameters of the N-cap off region observed 
by Suzaku with the emission measure of 0.12 cm“® pc (Tahara 
et al., 2015) which is at the Galactic longitude of 355.5 deg 
and the Galactic latitude of 35.8 deg, but we set the tempera¬ 
ture of 0.3 keV and the metallicity of 0.45 Z© for the bubbles. 
[O/Ee], and [Ne/Ee] of the bubbles are set to be zero, i.e. the 
solar abundance ratios. This situation roughly corresponds to 
the hadronic AW scenario. Under these assumptions, ASTRO- 
H/SXS can measure the metallicity of the bubbles as ZpB = 
0 . 45 lg' 2 iZ© and the abundance ratios as [O/Ee]= 0 . 00 lg ']^3 
and [Ne/Ee]=0.00to ii, where the errors represent 90% confi¬ 
dence level. If more metals are contained, metallicity and abun¬ 
dance ratios are more precisely constrained because of stronger 
line fluxes. Although precise determination of the metallicity 
is hard, we can determine abundance ratios precisely through 
the ASTRO-H observations. If ASTRO-H/SXS observe higher 
abundance ratios, it would strongly support the star forming ac¬ 
tivity scenarios as the origin of the bubbles. Moreover, precise 
determination of the abundance ratios will help us to distin¬ 
guish the gamma-ray emission process of the bubbles. 

4. Discussion and Conclusion 

In this paper, we showed that measurements of abundances 
in the bubbles will provide a unique clue to unveil their origin. 
The metal enrichment in the bubbles strongly depends on the 
bubbles formation scenarios and their emission mechanisms. 
It is still hard to determine the metallicities or abundances of 
the bubbles with current X-ray instruments. Eurther data or fu¬ 
ture missions are required. ASTRO-H/SXS can achieve a fac¬ 
tor of 10-100 times better energy resolution than SuzakufXlS 
do. Such high energy resolution will allow us to determine 
lines and their ratios. Based on the spectral simulation analysis, 
ASTRO-H/SXS will clearly detect line emissions. If high abun¬ 
dance ratios are obtained by ASTRO-H/SXS measurements, it 
will strongly support the star forming scenario as the origin of 
the bubbles. Moreover, precise measurement of the abundance 
ratios will enable us to investigate the gamma-ray emission 
process. Eurthermore, future X-ray mission Athena (Nandra 

^ Response files are taken from http: //astro-h. 

isas.jaxa.jp/researchers/sim/response. 

html. We adopt sxt-s_120210_ts02um_of_intallpxl.arf.gz 

for ARF, ah_sxs_7ev_basefilt_20090216.rmf.gz for RMF, and 

sxs_nxb_7ev_2011021 l_lGs.pha.gz for background files. 
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Energy (keV) 

Fig. 1. Simulated ASTRO-H/SXS spectrum of the Fermi bubbles with 
200 ks exposure. The data represents the expected performance by 
SXS, while the black, red, blue, and purple curve represents con¬ 
tributions from all components, the Fermi bubbles, the local hot 
bubble, and the cosmic X-ray background. We assume the spec¬ 
tral parameters of the N-cap off region observed by Suzaku with 
the emission measure of 0.12 cm“® pc (Tahara et al., 2015), but 
we set the temperature of 0.3 keV and the metallicity of 0.45 Zq 
for the bubbles. [0/Fe], and [Ne/Fe] of the bubbles are set to 
be zero, i.e. the solar abundance ratios. If more metals ex¬ 
ist in the bubbles, the stronger line emissions are expected. The 
position of the each line elements are indicated in the figure. 

et al., 2013) will have similar instrument but with higher en¬ 
ergy resolution and larger effective area. These future X-ray 
missions will enable us to understand the origin of the bubbles 
through the elemental abundances in the bubbles. 

For continuous injection models, we do not take into 
account the thermal conduction effect. As hot outflow gas 
exists behind the contact discontinuity, compressed gas can 
be heated up by the thermal conduction and flow behind 
the contact discontinuity. The abundance of the gas behind 
the contact discontinuity would be smaller than estimated. 
The thermal conduction time scale is given as tcond — 
10®(n/4 X 10“^ cm“^)(ZT/l-6 kpc)^(A:T/0.3 keV)“®'^^ yr 
(Kawasaki et al., 2002), where n is the gas density taken from 
Kataoka et al. (2013), It is the thermal conduction length 
assumed to be the thickness of the compressed region, and kT 
is gas temperature set to be 0.3 keV (Kataoka et al., 2013). 
Since the age of the bubble is expected to be in the order of 
10 Myr for the leptonic SF and leptonic AW scenarios, the 
results will not significantly change. However, in the case of 
the hadronic SF scenarios, the age would be comparable to the 
thermal conduction time scale. The actual abundance would 
be lower than that estimated in this paper. 

We assumed that the interior of the bubbles is described by 
a single temperature. In nearby starburst galaxies, observed 
X-ray emitting gas is composed of multi-temperature plasma 
(Strickland et al., 2002). Single temperature modelling may re¬ 
sult in erroneous abundance measurement. Here, the physical 
scale of the observed regions of the nearby starburst galaxies 
extends to ~ 3 kpc (Strickland et al., 2002), while that scale of 
the Field-of-View (FoV) of Suzaku/XlS and ASTRO-H/SXS at 
the GC is ~ 40 pc and ^ 7 pc, respectively. The expected fcond 
in the observable regions of the bubbles by Suzaku/XlS and 
ASTRO-H/SXS becomes much shorter than the age of the bub¬ 
bles. Thus, single temperature models work for the bubbles for 


pointing X-ray observations. Furthermore, the current X-ray 
spectra of the bubbles are well described by a single tempera¬ 
ture model (Kataoka et al., 2013; Tahara et al., 2015; Kataoka 
et al., 2015), although stacking analysis of the northern cap 
region indicates possible existence of another 0.7 keV plasma 
(Tahara et al., 2015). With ASTRO-H/SXS, we can observa- 
tionally distinguish another temperature component by com¬ 
paring the temperature based on single temperature spectral fit 
and that based on line ratios in each held. 

Non-thermal X-ray emission may underlie the thermal com¬ 
ponent as non-thermal emission is observed in radio and 
gamma-ray. Significant contribution of non-thermal emission 
may be crucial for deriving abundances. Kataoka et al. (2013) 
observationally constrained the non-thermal flux associated 
with the bubbles as < 9.3 x 10“® erg cm“^ sr“^ in the 
2-10 keV energy range, which is negligible comparing to the 
observed thermal flux. Theoretically, non-thermal X-ray flux 
of the bubbles is expected to be less than the observational 
upper limit through multi-wavelength spectral modelling (see 
e.g. Kataoka et al., 2013; Ackermann et al., 2014; Fujita et al., 
2014). 

We do not take into account the yields of Type la supernovae 
(SNe la) considering the uncertainties of the SNe la rate in 
the GC which is not observationally well constrained. SNe la 
are the thermonuclear explosions of accreting white dwarfs and 
produce Fe and little a-elements (e.g. Iwamoto et al., 1999). It 
is known that the cosmic SNe la rate is a factor of 3-10 lower 
than the cosmic core-collapse SNe rate (Horiuchi & Beacom, 
2010; Horiuchi et al., 2011). SN la explosion occurs not simul¬ 
taneously with star formation but delays. Delay time distribu¬ 
tion (DTD) of SNe la is represented by a power-law form (see 
e.g. Totani et al., 2008). By assuming a constant star forma¬ 
tion history (SFH) and a power-law DTD, the expected SNe la 
rate is ^ 0.01 per century which is roughly consistent with 
the estimate of 0.03 ± 0.02 per century (Schanne et al., 2007) 
based on the empirical relation between the rate and the stel¬ 
lar mass (Mannucci et al., 2005). The resultant iron abundance 
increases by 2% and 20% for the leptonic SF scenario (Lacki, 
2014) and the hadronic SF scenario (Crocker, 2012), respec¬ 
tively. We adopted the W7 model in Iwamoto et al. (1999) for 
the yields of SNe la and a power-law DTD following Yates 
et al. (2013). However, the SFRs in the nuclear bulge at ^ 30- 
70 Myr ago were about an order of magnitude lower than that 
at ~ 1 Myr ago (Matsunaga et al., 2011). Taking into account 
this SFH, the iron abundance does not change for the leptonic 
SF scenario, while it increases 2% for the hadronic SF sce¬ 
nario. Considering the uncertainties of future ASTRO-H/SXS 
measurements (see §. 3), the metal enrichment by SNe la in 
the bubbles would be negligible comparing to abundance mea¬ 
surement uncertainties. 

Abundances of stars and ISM in the GC region are assumed 
to be the solar. However, those abundances in the GC are still 
under debate. Various observations suggest that the GC metal¬ 
licity is at least in the range of ^GC $ 2 Zq (see the ap¬ 
pendix A of Crocker, 2012, for details), although their elemen¬ 
tal abundances are uncertain. If we assume 2Zq for ISM and 
stars in the GC, the resulting metallicity behind Red in the star 
formation scenarios increases by ~20% comparing to the case 
with solar abundance progenitor stars. We adopt the yields de- 
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scribed in Portinari et al. (1998) which give the yields for stars 
having up to 2.5^0, while the yields for stars having > Zq are 
not given in Nomoto et al. (2006). However, the effect of metal 
enrichment from HNs are not included in this comparison since 
those are not provided in Portinari et al. (1998). 

We thank the anonymous referee for useful comments and 
suggestions. We also thank to Noriko Yamasaki and Dmitry 
Malyshev for useful comments and discussions. Y.I. acknowl¬ 
edges support by the JAXA international top young fellowship. 
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